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We elucidate the relationship between Schrodinger-cat-like macroscopicity and geometric entan¬ 
glement and argue that these quantities are not interchangeable. While both properties are lost 
due to decoherence, we show that macroscopicity is rare in uniform and in so-called random physi¬ 
cal ensembles of pure quantum states, despite possibly large geometric entanglement. In contrast, 
permutation-symmetric pure states feature rather low geometric entanglement and strong and robust 
macroscopicity. 


I. INTRODUCTION 

Quantum entanglement entails two important conse¬ 
quences: On the one hand, it is the “characteristic trait 
of quantum mechanics” [1] that thoroughly thwarts our 
every day’s intuition, most bizarrely when applied to 
macroscopic objects, as illustrated by the famous para¬ 
dox of Schrddinger’s cat [2]. On the other hand, entangle¬ 
ment is the very ingredient that makes the simulation of 
quantum many-body systems extremely challenging: A 
separable system of N qubits requires only 2N parame¬ 
ters for its description, whereas an entangled state comes 
with ^ 2^ variables. This curse of dimensionality in the 
context of simulating quantum systems becomes a pow¬ 
erful resource when it comes to the speedup of quantum 
computers over classical architectures. Although the two 
aspects of entanglement are two sides of the same medal, 
they are quantified differently: Measures of macroscop¬ 
icity 0 reflect the degree to which a quantum state re¬ 
sembles a Schrodinger’s cat, the complexity of a quantum 
state is reproduced by the geometric measure of entan¬ 
glement [3], defined via the largest overlap to separable 
states. Macroscopicity can increase under local opera¬ 
tions combined with classical communication (LOCC) [3] 
in contrast to geometric entanglement, e.g. the modestly 
macroscopic cluster state can be converted via LOCC 
into a highly macroscopic GHZ-state [S]. 

In Nature, macroscopic entanglement does not occur 
P[7] outside artihcially tailored situations [5]. Its em¬ 
pirical absence contrasts with its immediate appearance 
in the formalism of quantum physics, which raises the 
question whether some unavoidable mechanism jeopar¬ 
dizes macroscopically entangled quantum states. The de¬ 
coherence programme [9] explains the emergence of clas¬ 
sical behavior in our everyday’s world and how macro¬ 
scopic quantum superpositions decohere on overwhelm¬ 
ingly short timescales The interaction between 

any quantum system and its surrounding environment 
destroys the coherence between macroscopically distinct 
alternatives - such as the dead and the living cat. Even 
if the environment were shielded off perfectly, however, it 
remains unlikely that macroscopic entanglement be ob¬ 
served in Nature, due to the unavoidable coarse-graining 


of any measurement [T4l - fT6] . In other words, there are 
powerful mechanisms that quickly deteriorate macro¬ 
scopic entanglement. 

But the decoherence programme does not make any 
statement on the likelihood that a macroscopic quantum 
superposition may form, only that, when it appears, it 
decoheres on a timescale so short that any attempt for ob¬ 
servation is vain. Would a hypothetical decoherence-free 
world then host cohorts of Schrodinger’s cats? In other 
words, how likely are macroscopically entangled quantum 
states before the onset of decoherence? 

In order to ease our intuition for these questions, 
let us propose a purely classical analogy: Consider 
the microcanonical ensemble of a gas of N particles 
in a box with total energy E, illustrated in Fig. 

In classical statistical mechanics, all microstates (spec¬ 
ified by the positions and momenta of the gas particles 
{xi,... ,xn]Pi, ... ,p]s[}) that are compatible with the 
total energy E are assigned the same probability m- 
How likely is it to find all particles in the left half of 
the box, as illustrated in Fig. [^b)? In the hrst place, 
if we prepared the particles in such a state, the system 
would relax to a homogenous distribution [Fig.j^a)] on a 
very short timescale - just like decoherence destroys any 
macroscopic quantum superposition. 
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FIG. 1: Two microstates of the microcanonical ensemble of 
2000 particles in a box. (a) Typical states feature a rather 
homogeneous distribution of particles, while (b) states with 
large inhomogeneities are artificial and rare. We argue in 
this article that the rarity of macroscopically entangled pure 
quantum states can be understood in close analogy: Macro¬ 
scopically entangled states play the role of states with large 
inhomogeneities (b), and therefore form only extremely sel¬ 
dom. 
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This explanation for the absence of the strongly 
inhomogeneous situation is, however, not the com¬ 
monly adopted one: The spontaneous occurrence of the 
macrostate sketched in Fig. 0b) is by itself extremely 
unlikely, and the relaxation described above is in fact 
not required to explain its rarity. There are overwhelm¬ 
ingly more microstates that correspond to a homogeneous 
distribution (a) than for a distribution with all parti¬ 
cles on one side (b) [T^. We can therefore safely ne¬ 
glect the inhomogeneous macrostate and focus on the 
macrostate with a homogeneous distribution of particles 
- which is the very basis of statistical mechanics m- 
In other words, even though there is a powerful mech¬ 
anism to restore the homogeneity of the gas particles, 
this mechanism is not required to explain the absence of 
inhomogenous distributions: These macrostates are suf¬ 
ficiently unlikely to occur to be safely neglected. Quite 
similarly, as we will argue below, despite decoherence be¬ 
ing a powerful mechanism to explain the disappearance 
of macroscopic quantum superpositions uni n], these 
states are a priori extremely rare in many ensembles of 
pure states. 

While mixed states already include the effect of de¬ 
coherence and feature little macroscopicity [3], we fo¬ 
cus here on pure states. The typicality of macroscopic 
entanglement then depends on the actual choice of the 
pure-state ensemble. We find analytical and numeri¬ 
cal evidence that macroscopic superpositions are untyp¬ 
ical among random pure states in different ensembles. 
This finding may appear paradoxical at first sight, since 
macroscopicity leads to multipartite entanglement [TS] 
and entanglement, as quantified by the geometric mea¬ 
sure in [19] (tantamount to a large distance to separa¬ 
ble states), is common in random pure states [IHl I3T] . 
We resolve this ostensible paradox by showing that geo¬ 
metric entanglement is actually adverse for macroscopic¬ 
ity: Random states are very entangled, and hence non- 
macroscopic. 

We review measures of macroscopic and geometric en¬ 
tanglement in Section [ll| where we also present some 
technical results. Since both quantities are dehned via 
a maximization procedure, their relationship is intricate. 


We elucidate their connection qualitatively in Section III 


in order to gain a good understanding of the statistics 
of macroscopic entanglement evaluated in different pure- 
state ensembles in Section |IV| We conclude in Section jV] 
where we propose an extension of our study to other en¬ 
sembles and sketch its consequences for the preparation 
of macroscopic states in the experiment. 


II. QUANTIFYING MACROSCOPIC AND 
GEOMETRIC ENTANGLEMENT 

In order to address our central question - is macroscop¬ 
icity typical in ensembles of pure quantum states? - we 
need to establish a quantitative measure for macroscop¬ 
icity. Our system of interest is a collection of N qubits, 


which, for the purpose of illustration, we treat as spin 
1/2-particles. That is, we focus on pure quantum states 
living in the Hilbert space H = . While no consen¬ 

sus exists on how to rigorously quantify macroscopicity 
for mixed states [3l|22H27], this debate is not crucial in 
our context, since we focus on pure states, for which most 
measures agree. We will adapt a well-established measure 
for macroscopic entanglement [ 3121 ] and propose a sensi¬ 
ble way for its normalization. To set the context, we are 
interested in the typicality of macroscopic entanglement 
within quantum theory; a general benchmark of macro¬ 
scopic quantum superpositions will require concepts that 
draw beyond this realm [29l [30] , and may be tailored for 
specific applications m- 


A. Measure of mascroscopicity 

Macroscopicity manifests itself in disproportionally 
large fluctuations of some additive multi-particle observ¬ 
able [3132], i.e. of some operator of the form 

N 

s{di,...dN) = '^dj-dj, ( 1 ) 

where dj is the vector of three Pauli matrices that act on 
the jth qubit and the local orientation of the measure¬ 
ment operator dj is normalized, 

|a,p = l. (2) 

The operator S describes the total spin of the system 
with respect to locally adjusted spin-directions, defined 
by dj = Given a pure quantum state |'I'), 

the maximally obtainable variance of this additive oper¬ 
ator, 

{AS^{d^,...dN)) = (vI/|^2|vI/)-((vI-|5|M'))2, (3) 

defines the unnormalized macroscopicity M [3, 

M = max5i,...^5„(A,S^(di,.. .a^v))- (4) 

Quite naturally, we find 

N<M<N^. (5) 

The lower bound is saturated, e.g., for separable states, 
for which we maximize fluctuations by choosing the mea¬ 
surement directions dj to be unbiased with respect to the 
local spin directions. The upper bound is reached, e.g., 
for the Greenberger-Horne-Zeilinger (GHZ) state [33] 

|GHZ)^r = ^(|0)®^ + |l)®^), (6) 

which comes closest to modelling a coherent superposi¬ 
tion of a living and a dead cat: The two super-imposed al¬ 
ternatives (all spins pointing up, all spins pointing down) 
are maximally different, leading to maximal fluctuations. 
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In order to compare the macroscopicity of systems of 
different sizes iV, we normalize Ai: 


/jCi(|«))-JV 

such that 0 < AI < 1 for all TV, where the upper (lower) 
bound is saturated for M = (N). 


B. Additivity and basic properties 

The unnormalized macroscopicity A4 is additive in the 
sense that any product state |'h) (8) |$) yields 

AI(|^)®|$))=Ai(|^))+7W(|$)), (8) 

because the separability of Idt) (g) |$) excludes additional 
fluctuations by choosing a direction of spins other than 
the optimal ones for Idt) and |$). 

Any family of states for which the size of non-separable 
components does not scale with the system size has van¬ 
ishing normalized macroscopicity in the limit of many 
particles. For example, a tensor product of N/2 Bell- 


states is not macroscopic: 




= 2N 

(9) 

Ai 

1 N—>-oo ^ 

= , 0 

(10) 


C. Relation to index p 

Our measure of macroscopicity can be directly related 
to the index p [33] , the exponent that defines the scaling 
of A4 with N, 

M(xNP. (11) 

The scaling properties of macroscopic entanglement can 
be investigated for state families, i.e. “prescriptions that 
assign to any system size a quantum state |'h)” [32] • Any 
state family for which A4 > S > 0 with S independent of 
N can be considered macroscopic: In the limit N —> oo, 
A1 is then related to the index p as follows: 

M = 0 ^ p = l, (12) 

7W>(5>0 p = 2. (13) 

In other words, for p = 2, we have macroscopically large 
fluctuations that do not vanish in the limit of many parti¬ 
cles; the index p, however, does not yield any information 
about the actual fraction of particles participating in a 
macroscopic superposition. The normalized macroscop¬ 
icity A4 is more fine-grained and offers such insight: For 
the family of states 


(14) 


we have 


TWd^-)) = Nf + N2, 


and, thus. 


M(|vE-)) 


Afi,JV2>l 


Ni 

Ni + N 2 


(15) 


(16) 


In this case A4 reflects the fraction of particles in the sys¬ 
tem that take part in a quantum superposition of macro¬ 
scopically distinct alternatives. On the other hand, for 
large N, we also have NAi « vCm (with an additive 
error of the order ~ VN), which quantifies the absolute 
size of the macroscopic superposition. In the following, 
we will therefore focus on the macroscopicity Ai, as de¬ 
fined in Eq. Q. 


D. Evaluation of macroscopicity 

1. Generic states 

The definition of the unnormalized macroscopicity Ai 
in Eq. Q entails an optimization problem over 2N vari¬ 
ables, which evidently complicates its evaluation for large 
systems. As an alternative to multivariable optimization, 
it was proposed [33] to evaluate the unnormalized macro- 
copicity A4 using the variance-covariance-matrix 

V^k,pj = Ivh), (17) 

where AaJ = aj — (dp'), 7 = x,y,z. The resulting 
3A X 3A-matrix V then stores the fluctuations of all 
observables that are sums of Pauli matrices, and we have 

m 


N 

(A5^(ai,...,aAr)) = ^ ^ aJVoy/jfeof. (18) 

j^k—1 'y,(3—x,y,z 

The expectation value (A5'^(ai,... ,aAr)) is maximized 
by choosing the dj as the Eigenvector vi of V correspond¬ 
ing to the largest Eigenvalue Ai of E, such that 


Yj = VN I Wi,3Q_i)+2 
V Vl,3U-l)+3 


(19) 


However, the aj chosen this way are only constrained by 


N 


J2\di\^ = N, 


( 20 ) 




which is a much weaker constraint than our Eq. Q: 
Instead of N unit-normalized Bloch-vectors dj, only 
the sum of the norm of all Bloch-vectors is fixed in 


Eq. (20). Colloquially speaking, the dj chosen accord¬ 
ing to Eq. (19) [33] allow us to weight the importance 
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of individual spins in the system differently, and give 
those featuring large fluctuations a larger impact. We 
can therefore only state that 

^(|«'))<iVAi. (21) 

Admittedly, the scaling of Ai with the system size N for 
a given family of states is inherited by A4 evaluated via 
Eq. such that coarse-grained quantities such as the 
index p can be evaluated using Ai. In particular, the 
Eigenvalues of V only take the values 1 and 0 for sepa¬ 
rable states, which are clearly - and not surprisingly - 
not macroscopic at all. For a quantitative understanding 
of macroscopicity, however, the optimization inherent to 
0 is crucial: For example, consider the state 

l^c) = |^+) 0 |0)®^-", (22) 

with N > where the first two qubits are in a maxi¬ 
mally entangled Bell-state, but remain completely sep¬ 
arable from the rest of the system. The unnormalized 
macroscopicity of |$c) fulfils 

M=N + 2<2N = XiN, (23) 

i.e. the largest Eigenvalue Ai is related to the Eigenvector 
ui, with 


2. Symmetric states 


In the case of permutation-symmetric states, V as¬ 
sumes a structure with repeated 3 x 3-blocks, and the 
Eigenvectors Vj reflect this symmetry. As a consequence, 
the optimal local spin orientations all coincide, a/ = dk 
for all k,l, and, consequently, Pj = dj. The optimiza¬ 
tion inherent to Eq. Q becomes unnecessary, since the 
lower bound (281 and the upper bound (21) on the unnor¬ 
malized macroscopicity coincide. We can therefore safely 
adopt the method introduced in Ref. [35] to compute A4. 

Specifically, the variance-covariance-matrix V then 
consists of two different 3 x 3-blocks 


(29) 

(30) 


which contain all variances and co-variances, respectively, 
and assumes the structure 


V = 


/ A B ... B\ 
B A ... B 


\B B ... A J 


(31) 


(X\ — 

rN 

Vy(1.0,0) 

(24) 

Oi2 = 

Vy(1.0,0) 

(25) 

dk>d. = 

(0,0,0) 

(26) 


which is not compatible with Eq. ([^ and entirely ignores 
the separable qubits 3,... ,N, while the two entangled 
qubits are over-weighted. This example being admittedly 
artificial, we have nevertheless experienced a substantial 
difference between the exact calculation and the value 
extracted via the VCM for the non-symmetric ensembles 
of random states considered below in Section 

O n th e other hand, given a set of orientations obtained 
via (19) and \djp > 0 for all j, we can normalize the dj 
to find a candidate spin orientation that promises to yield 
a large variance. 


By writing V = 1 ^ A + M ^ B, one can show that 


the largest Eigenvalue Ai of the above block-matrix V 
coincides with the largest Eigenvalue of the 3 x 3-matrix 

V,y^ = A+{N-l)B, 

(32) 

such that 


Ad(|^,y^)) = iVAl. 

(33) 


This greatly facilitates the computation of the macro¬ 
scopicity for permutation-symmetric states. To obtain 
the matrices A and B, we use the efficient techniques for 
the computation of reduced density matrices of symmet¬ 
ric states presented in Ref. |5B] . 


E. Geometric measure of entanglement 


/ 3 .= 




(27) 


The resulting value (A^^(/3i,...,/3jv)) then provides a 
lower bound on the actual unnormalized macroscopicity, 
since we are not guaranteed that the choice of local spin 
orientations given by (|27|) is the optimal one: 


We will relate macroscopicity to the geometric measure 
of entanglement HISzl, which is defined via the maxi¬ 
mal overlap ij of |4>) with any separable state l^sep) = 
\<Pij4'2, ■ ■ ■ , Pn), 


Eg(|^')) =-log 2 ?? =-logs sup |($sep|4')r (34) 

I'f.ep) 




(28) 


In other words, even though the precise value of Ai re¬ 
quires a numerical optimization, computationally inex¬ 
pensive lower and upper bounds [Eqs. (28) and (21), re¬ 
spectively] to this quantity can be established straight¬ 
forwardly. 


The geometric measure of entanglement of N qubits nat¬ 
urally vanishes for separable states, and is bounded from 
above by N — 1. High geometric entanglement is tanta¬ 
mount to a large generalized Schmidt measure |38j , which 
reflects the number of separable terms required to express 
the state; the geometric measure of entanglement, thus, 
quantifies the complexity of a quantum state. 
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Since we will face the statistics of geometric entangle¬ 
ment in Section IV we discuss its evaluation in practice. 
In general, the computation of the geometric measure of 
entanglement requires an optimization over the 2N free 
parameters xi,... ,xn and yi,... ,yN that define the sep¬ 
arable state 


l^’sep) = (cosa;j|0) -|- sinxj|l)) , (35) 

entailing signihcant computational expenses. Alterna¬ 

tively, candidate solutions for the closest separable state 
can be computed via the probabilistic iterative algorithm 
presented in Ref. [TU] . 

For permutation-symmetric states, the evaluation of 
the geometric measure is facilitated considerably. A 
permutation-symmetric state of N qubits can be writ¬ 
ten in the Majorana-representation, 

^ E (36) 

o'G.S N 

where 

Af = N\ perm((ej|efe)), (37) 

is a normalization constant, and perm((ej|efc)) is the per¬ 
manent [39] of the N X N Gram-matrix that contains 
all mutual scalar products (ej|efc). The permanent is, 
in general, a function that is exponentially hard in the 
matrix size N, but for a Gram-matrix with only two non¬ 
vanishing singular values, as given here by construction, 
an efficient evaluation is possible [39] . For this purpose, 
the representation of symmetric states in the Dicke-basis 
is valuable. 


N 

(38) 

j=0 

The Dicke-states are defined as 

/ AJ\ 

l^!^^)=(.) E ( 39 ) 


where the summation includes all possibilities to dis¬ 
tribute j particles in |I) and N — j particles in |0) among 
the N modes. The quantitative relationship between the 
expansion coefficients in the Dicke-basis Cj and the N 
states |efc) that define the Majorana representation (36) 
is presented in Ref. HO). 

Since the closest separable state to a symmetric state 
is itself symmetric |41] . the optimization problem over 
2N variables implicit in ( |34| ) reduces to a merely two- 
dimensional setting. Given the Majorana-representation 
|ei),..., I Eat), we need to find the single-qubit state {(j)) 
that maximizes 


1 ^ 

( 40 ) 


Besides standard numerical optimization strategies, we 
can adapt the iterative algorithm of Ref. [19] to symmet¬ 
ric states: For that purpose, we choose a random single¬ 
particle state |(/)o). We then iteratively generate states 
\(j)k) as follows: 


<//c+l) 

AL — 

II 

(41) 

<//c+l) 


(42) 


{^k+i\^k+i) 


The results of Ref. [19] can be translated to a wide extent 
to the current setting with symmetric states, and 
becomes a good candidate for the closest separable state 
for large fc, although the algorithm is prone to return a 
local instead of a global maximum. 


III. RELATIONSHIP BETWEEN 
MACROSCOPICITY AND GEOMETRIC 
ENTANGLEMENT 


Having established the two pertinent main character¬ 
istics of a quantum many-body state - its macro- 
scopicity Ai [Eq. 0] and its geometric entanglement Eq 
[Eq. (@], we can now explore the relationship between 
these two quantities. 

The behavior of entanglement measures and macro- 
scopicity was analysed in Ref. [T5] , where the author hnds 
that “a state which includes superposition of macroscop- 
ically distinct states also has large multipartite entangle¬ 
ment in terms of the distance-like measures of entangle¬ 
ment” . Here, we will confirm that geometric entangle¬ 
ment is indeed necessary for non-vanishing macroscopic- 
ity; however, we will also show that the general relation¬ 
ship between macroscopicity and geometric entanglement 
is rather involved. In particular, non-vanishing geometric 
entanglement is not sufficient for macroscopicity: There 
are entangled states with strictly vanishing macroscopic¬ 
ity. On the other hand, very large geometric entangle¬ 
ment implies small macroscopicity, i.e. the maximal value 
of macroscopicity is reached for finite geometric entangle¬ 
ment. As a consequence, the two quantities should not 
be used synonymously. 


A. Close-to-separable states 


We first focus on states |4>) for which the largest 
squared overlap with a separable state ry is larger than 
or equal to 1/2, i.e. the geometric measure of entangle¬ 
ment Acd'k)) [Eq. (34)] is smaller than or equal to unity. 


To this end, we explore the transition between the sep¬ 
arable state jo, ..., 0), which carries neither macroscop- 
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icity nor entanglement, to the maximally macroscopic 
GHZ-state (|^, described by 


|S(0,c)) 


cos6\0)'^^ + sin0(cose|O) + sine|l))®'^ 
+ cos^ e sin(20) 


with 0 < e < 7r/2 and 0 < 0 < 7r/4. The state |S(6>, e)) 
is a superposition of two separable components in which 
all qubits populate the very same states, with weights 
depending on 9. The distinguishability of the two al¬ 
ternatives is defined by e. The parametrization in e for 
fixed 6 = tt/4 was introduced in Ref. [12] and explored in 
Refs. 131421 . Since the state is permutation-symmetric, 
we can use the methods of Section HID 21 to evaluate its 
macroscopicity. 

Several limiting cases are reached for particular values 
of the parameters 9, e: For {9, e) = (7r/4,7r/2), we recover 
the GHZ-state ([^; for 0 = 0 or e = 0, we deal with a 
separable state. For e > n/2, the destructive interference 
between the two amplitudes associated to the component 
|0)®'^ can lead to geometric entanglement larger than 
unity. In particular, we obtain the IT-state 



Geometric entanglement Eg{\=l{6, e))) 


FIG. 2: Macroscopicity as a function of geometric entangle¬ 
ment for the family of states |H(0, e)). The two extremal cases 
are given by e = 7r/2 (black dashed) and 0 = 7r/4 (solid red), 
which take turns as the upper bound on macroscopicity for 
given geometric entanglement. 


the pertinent variance-covariance matrix I4ym 
then becomes 


in Eq. (32) 


|IF) = ID^) = ^ (|1,0,..., 0) + |0,1,0,..., 0) 

in the limit e —>■ tt, 0 —?> 7r/4, for odd N. 

We parametrize the closest separable state as 


)(44) 


/I 0 0 \ 

= 0 1 0 

\ 0 0 N 3111^29 j 


which yields the macroscopicity. 


<l>sep(a)) = (cosa|0)-I-sina|I)) 


®N 


(45) 


the overlap with |i:.(0, e)) becomes 

|($.ep(a)|S(0,e))|2 = 


2 (cos0cos^ a-b sin0cos^(e — a))^ 


M = 


I max(l, fVsin^ 29) — 1 


1 -b cos^ esin(20) 


{,46) 


1 

The geometric entanglement takes the value 
Eg = -log2(cos^6i). 


(47) 


(48) 


(49) 


which needs to be maximized with respect to a to obtain 
the geometric measure of entanglement. For large the 
overlap is maximized for a = 0 (since 9 < 7r/4); for finite 
N^ the maximum can conveniently be found numerically, 
since the overlap (461 does not oscillate fast as a function 
of a. 

We show the behavior of geometric and macroscopic 
entanglement in Fig. [^ for different numbers of qubits 
N. Although macroscopicity increases with geometric en¬ 
tanglement as a general tendency, the relationship is am¬ 
biguous, especially for large numbers of qubits N. Based 
on extensive numerical evidence, we conjecture that the 
general maximum macroscopicity for a given value of geo¬ 
metric entanglement is attained by the value obtained for 
|S(0 = 7r/4, e)) or |5(0,e = 7r/2)), in the range Eq < 1. 


1. Geometric entanglement without macroscopicity 

We first consider the family of states parametrized by 
e = 7r/2,0 < 9 < tt/4 (black dashed lines in Fig. [^, 


As a consequence, for 0 < 0 < | arcsin yJl/N, the geo¬ 
metric entanglement remains finite, yet the macroscop¬ 
icity vanishes. In other words, there are states that are 
entangled, but the fluctuations in any additive observable 
do not surpass those that can be achieved for separable 
states. This explains the step-like behavior observed in 
Fig.H best visible for small N. On the other hand, 
even though geometric entanglement is not sufficient for 
macroscopicity, it is necessary: For a separable state, all 
non-vanishing Eigenvalues of the variance-covariance ma¬ 
trix 0 are unity. 


2. Close-to-unity geometric entanglement and small 
macroscopicity 

The most macroscopic state, the GHZ state pos¬ 
sesses geometric entanglement £’g(|4>ghz)) = 1- An¬ 
ticipating Section Ell large macroscopicity naturally 
comes with geometric entanglement Eq ~ I, in general. 
The criterion Eq « 1 is, however, not sufficient to en¬ 
sure large macroscopicity: Consider the family of states 
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6 = 7r/4 



0.8 


o.e 


0.4 



0.0 0.1 0.2 0.3 0.4 0.5 

e/TT 


FIG. 3: Decomposition of |H(0,£)) for N = 2A into Dicke- 
state-components [Eq. (39l], for the two extremal families of 
states characterized by e = 7r/2 (upper panels) and 0 = 7r/4 
(lower panels). For e = 7r/2, only the very first and very 
last Dicke-components are populated, while a binomial dis¬ 
tribution of components slowly shifts to the highest Dicke- 
component for 9 = 7r/4. Consequently, the behavior in the 
Eg — A4-plane [Fig.[^ is very different for the two families of 
states. In the lower right panel, the black dashed line shows 
the upper bound to Eq given in Eq. (511. 


parametrized by 0 = 7 r/ 4,0 < e < 7 r/ 2 . By evaluating 
the largest Eigenvalue of the matrix I4yni, we obtain the 
macroscopicity 


M 


■ 2 
sm e 


1 -I- cos 


N , 


(50) 


The maximal overlap to separable states is bounded from 
below by the overlap with the separable test-state (45) 
setting a = 0 ; the geometric measure of entanglement 
therefore fulfils 


is separable. For decreasing e « 7 r/ 2 , we continuously 
loose macroscopicity, because the average directions of 
the spins become similar and the two superimposed al¬ 
ternatives less and less macroscopically distinct. How¬ 
ever, the closest separable state remains the Dicke-state 
|I?®) for a wide range of e < 7 r/ 4 , i.e. the geometric 
measure of entanglement remains close to unity. In con¬ 
trast, for the family parametrized by e = 7r/2, the loss 
of geometric entanglement is directly accompanied by a 
loss of macroscopicity (upper panel of Fig.[^. 


B. Far-from-separable states 

Let us now move into the domain of strongly geomet¬ 
rically entangled states and assume that we are given a 
maximal overlap with separable states i] < 1/2, i.e. a ge¬ 
ometric measure of entanglement Eq = — log 2 ry > 1. We 
construct a state that maximizes the macroscopic¬ 
ity under this constraint. Since local rotations do not 
affect the geometric measure of entanglement, we can as¬ 
sume that the optimal value of the spin-orientations all 
point into the z-direction (oj = (0, 0,1)), and we expand 
Ifh) in Eigenstates of the total spin operator S, which 
has Eigenvalues —N, —N + 2,..., N. Each Eigenvalue 
is ((Ar-i-S)/ 2 )"^°^®^ degenerate, the state can therefore be 
written as 

((S + N)/2) 

E E cs.a|5,A), (52) 

S=-N-N+2,...N A=1 

where A labels the degenerate states, and we choose the 
l^, A) to be separable. The expectation value of powers 
of the collective spin becomes 

((S + N)/2) 

= E E |cs,Ap.(53) 

S=-N,-N+2,...N A=1 

Under the constraint that the maximal overlap to any 
separable state be fixed to rj, 

|cs,Ap < ??, (54) 


f 1 _1_ 

EG{\E{e = 7r/4,e)}) < - log^ (51) 


That is, for a wide range of e < 7 r/ 2 , we retain a geometric 
measure of entanglement of around unity, but quickly 
loose macroscopicity. In Fig. the red lines quickly dive 
into low values of macroscopicity, while remaining close 
to Eq = 1, a trend that becomes more and more clear 
for larger values of N. This stands in stark contrast to 
the black dashed lines, which retain macroscopicity for 
decaying geometric entanglement. 

This behavior can be understood intuitively via the 
decomposition of the state into Dicke-states [Eq. (38)], 
shown in Fig. The largest overlap with any sepa¬ 
rable state is at least as large as the coefficient in the 
Dicke-state expansion related to since the latter 


we maximize the variance ( 1 ^ by setting cn^i = C-n,i = 
1^/77 and subsequently distributing the probability ampli¬ 
tude among the remaining coefficients, i.e. we 

set as many pairs of coefficients to cs,k = C-s,k = \/v 
as possible, proceeding from large to small total spins S. 
That is, we maximize the contribution to the expectation 
value of S^, while the expectation value of S remains 0. 
The last pair cs,k = C-s,k is set to accommodate the re¬ 
maining amplitude, typically smaller than yd). Formally, 
the state reads 

Ill^T^^niax) ~ ^ ^ 

S=-N,-N+2,...,-nwd(N,2) 

Un-S)/^) 

Y, Vv k) + I - S, k)) , (55) 

k=l 
















A. Random physical states 


where the sum only runs over so many terms such that 
the state is normalized to unity, one term may possibly 
be weighted by a factor smaller than 


The obtained bound is shown in Fig. I^for different val¬ 
ues of N, as a function of geometric entanglement. The 
maximally achievable macroscopicity grows as a function 
of N for a fixed overlap 77 , but, for a fixed number of 
qubits N, a small overlap 77 , equivalent to large geomet¬ 
ric entanglement, causes a reduced macroscopicity. The 
expansion into separable states Eq. (521 is, however, not 
necessarily the optimal generalized Schmidt decomposi¬ 
tion [^, i.e. it is often possible to find a separable state 
with overlap larger than 77 . We can therefore not expect 
the bounds to be tight. 


We can repeat the argument for symmetric states, for 
which we impose that the amplitude of each Dicke-state 
component is constrained by 77 , leading to a state 

of the form 




0 = 


E Vv 

k=0...[N/2\ 


\D 


(fc)\ 

N / 




Due to the strict symmetry constraint on the state, we 
obtain smaller maximal values of macroscopicity for given 
geometric entanglement (dashed lines in Fig. @) than 
for general states. 



FIG. 4: Maximal value of normalized macroscopicity, evalu¬ 
ated using the argument of Section mE for general states 
(solid lines) and symmetric states (dashed lines), for N = 
3, 6,10, 20, 30 (blue, red, black, green and orange, respec¬ 
tively) . The horizontal dotted line in dicates the limiting value 
l/%/3 for symmetric states (Section IVCl. For A = 3, the 
bounds for general and for symmetric states coincide. 


IV. STATISTICS OF MACROSCOPIC AND 
GEOMETRIC ENTANGLEMENT 


Having established the relationship between geometric 
and macroscopic entanglement, we proceed to numerical 
investigations of pure states in different ensembles. 


1. State generation 

As a first ensemble of pure states, we consider random 
physical states, introduced in Ref. [43l |44j and sketched 
in Fig. I^a). To generate a random physical state [dlfc), 
we assume that the qubits are aligned in spin-chain con¬ 
figuration and that only pairwise interactions take place. 
We apply k times a random two-particle unitary onto a 
randomly chosen pair of two neighboring qubits. That is, 
for fc < A — 1 , the state remains at least 1 -separable (the 
first qubit in the chain has never directly or indirectly in¬ 
teracted with the last one), while we obtain Haar-random 
states in the limit k ^ N, which we will discuss sepa¬ 
rately in Section pV B| below. To obtain a better intuition 
for this ensemble, we compare random physical states to 
random linear chains |$fc)j which are generated by ap¬ 
plying a binary unitary between the first k < N —1 pairs 
of qubits, starting from a separable state [Fig. [^b)]. As 
the authors of [43l 04] argue, the ensemble of random 
physical states can be regarded as typical for physical 
systems that obey some locality structure. A variant 
of random physical states for which the binary interac¬ 
tions are chosen to be not necessarily between adjacent 
neighbors but between any two randomly chosen qubits 
does not exhibit qualitative differences to the locality¬ 
preserving model here. 


S. Numerical results 

Geometric and macroscopic entanglement for random 
physical states are shown in Fig. [^c,e) as a function 
of the number of applied binary gates k, which can 
be confronted to the behavior of random linear chains 
[Fig. [^d,f)]. Between fc « A and k « A^, we observe a 
steep increase in the geometric entanglement of random 
physical states: For this range of numbers of binary inter¬ 
actions k, the state typically becomes fully inseparable. 
For k K. A^, we observe a saturation of both macro¬ 
scopicity and geometric entanglement. While geometric 
entanglement increases monotonically with the number 
of applied gates, macroscopicity develops a peaked struc¬ 
ture for A > 6 . Consequently, the trajectories of random 
physical states as a function of k in the (Al, Acj-plane 
[Fig.§ proceed from low geometric and macroscopic en¬ 
tanglement (fc = 1 ) over a maximum to the asymptotic 
value with large geometric and low macroscopic entangle¬ 
ment. The maximum value of macroscopicity is reached 
for A < fc < A^: In this range, the states can very proba¬ 
bly not be decomposed into separable components, while 
it remains moderately complex by construction - these 
are the very requirements for high macroscopicity. Ran¬ 
dom linear chains feature a linear increase of geometric 
entanglement with the number of applied gates, for which 
the curves for all particle numbers coincide [Fig. [^d)], 
and a monotonic increase of macroscopicity for fc, peak- 
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ing at lower and lower values as we increase the number 
of particles [Fig. [^f)]. 


= 0 ^ ^ ^ ^ 

k = if^ 






Applied binary gates log^(/c) 





k = 0 




k = 1 


k = 2 





0 5 10 15 

Number of applied gates k/(N-l) 


(f) 


^ 0 ., 


V=4 


■fl 




j J.yiH* 




Z 0.0 0.2 0.4 0.6 0.8 1.0 

Number of applied gates k/(N-l) 


FIG. 5: (a) Random physical states. We apply k times a 
random unitary between two randomly chosen adjacent sites 
(closed boundary conditions). For k ^ N, we have a fully 
connected system with high probability, i.e. the state is typ¬ 
ically 0-separable, for fe —>■ oo, we reach the limit of Haar- 
uniform states, (b) Random linear chain. We apply random 
unitary binary gates between the first fc < — 1 pairs of adja¬ 

cent qubits. For fc = A^ — 1, we have a fully connected system. 
Geometric entanglement (c,d) and normalized macroscopicity 
(e,f) for random physical states and random linear chains, for 
k random two-body unitaries, respectively. In (d), the geo¬ 
metric entanglement coincides for all numbers of qubits N. To 
compare different system sizes, we plot the entanglement and 
the normalized macroscopicity as a function of the A'’-base 
logarithm of k in (c,e), and as a function of the normalized 
number of applied gates in (d,f). Sample size is 200, error 
bars show one standard deviation. 

Macroscopicity as a function of the particle number N 
is plotted in Fig.[^ The maximum value in random phys¬ 
ical states decreases with increasing N (red diamonds), 
albeit slower than the saturated value of the macroscop¬ 
icity (fc = N^, black circles). The latter remains slightly 
lower than the macroscopicity reached for a saturated 
random linear chain (i.e. after k = N — 1 binary gates) 
- choosing the interacting qubits randomly is disadvan¬ 
tageous for large macroscopicity, which gives an advan¬ 
tage to saturated linear chains. Complexity is adverse to 



FIG. 6: Average trajectories of random physical states 
in the Sg — Al-plane for = 4, 6,13, 20. The solid lines start 
at fe = 1 and proceed to fc = (blue discs). 




FIG. 7: (a) Normalized macroscopicity M after one random 

gate (fc = 1, green squares), after k = random gates 
(black discs), maximal attained macroscopicity for random 
physical states (red diamonds) and normalized macroscopicity 
for a binary chain with k = N — 1 random interactions (blue 
triangles), as a function of the number of qubits N. (b) Scaled 
macroscopicity NM ~ M for the same ensembles of states 
as in (a). Error bars show one standard deviation. 


macroscopicity: For N > 15, one randomly chosen binary 
interaction onto two qubits in an initially separable state 
[k = 1, green squares) results in a larger macroscopicity 
than the limiting case k —>■ N^. 

The different types of decay raise the question whether, 
albeit the fraction of particles participating in macro¬ 
scopic superpositions decreases, the absolute number 
may in fact be constant or increase. We plot the absolute 


size of the macroscopic component NA4 « - the ap¬ 

proximation is justified for A >> 1 which resolves the 
qualitative differences between the ensembles: The abso¬ 
lute number of particles participating in a GHZ-like state 
remains constant for states into which exactly one ran¬ 
dom gate has been applied (green squares), it decreases 
for k = (black circles), but it increases with N for 
the maximally achieved value in random physical states 
and for saturated random linear chains {k = N — 1). 

In conclusion, starting from a separable state and ap¬ 
plying random binary gates, we first explore the region in 
which geometric and macroscopic entanglement are syn¬ 
onymous (Section III A I, such that both quantities ini- 
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tially grow with k. When the spin-chain is fully insepara¬ 
ble, additional interactions contribute to larger geometric 
entanglement, but simultaneously destroy its mac roscop - 
icity ensuring that the latter decreases (Section IIIB). 
Even thought the absolute size of the macroscopic com¬ 
ponent increases with N [Fig. Bb)], the fraction of par¬ 
ticles participating in a macroscopic superposition does 
not. To come back to our classical analogy [Fig. 0 , just 
like there are no concerted forces that spontaneously push 
all gas particles to one side of the box, random evolutions 
are unlikely to force all spins into a macroscopic super¬ 
position. 


B. Haar-random states 

1. State generation 

In the limit fc —> oo, random physical states converge 
to Haar-random states, i.e. the ensemble of pure quan¬ 
tum states that are uniformly distributed on the unit 
sphere in Hilbert space [45) . Instead of applying many 
binary gates, one can construct Haar-random states by 
randomly generating the real and imaginary part of each 
state coefficient following a zero-mean unit- 

variance normal distribution, the resulting unnormalized 
vector c is then normalized in a second step. This pro¬ 
cedure yields a “chaotic” ensemble [46] that remains in¬ 
variant under local basis-rotations and re-partitioning of 
the Hilbert-space into subsystems m- Random states 
also result from the application of a Haar-random uni¬ 
tary matrix on any constant pure state [48] . 


2. Macroscopicity is rare in Haar-random states 

Random states feature the concentration of measure 
phenomenon, i.e. most states on the high-dimensional 
Bloch-sphere lie close to the equator |49|. Given a 
Lipschitz-continuous function /([dt)), the function values 
remain close to the average value (/) for the vast major¬ 
ity of states, reflected by the probability for a deviation 
larger than e [50] . 

P[\fm)-{f)\>e]<4e-^^, (57) 

where rj is the Lipschitz constant. Using trial func¬ 
tions, we find that the macroscopicity defined in Eq. Q 
is Lipschitz-continuous, while the geometric measure 
of entanglement inherits Lipschitz-continuity from the 
distance-like measure it is based on |33|. Hence, most 
Haar-random states are very similar, both when char¬ 
acterized by their geometric entanglement and by their 
macroscopicity. 

Using random matrix theory m , one can estimate 
the typical magnitudes of the elements of the variance- 
covariance-matrix ( [T7| ) [46] . In the limit N —>■ oo, the 
VCM approaches the unit matrix and, as the largest 


Eigenvalues converge to unity, by the upper bound (211, 
the normalized macroscopicity vanishes. This result 
also follows from the following complementary argument: 
Random states that are chosen according to the Haar 
measure possess large geometric entanglement 12011111: 
With probability greater than 1 — , we have for 

TV > 11 1201 


£;G(|^random))>A^-2l0g2TV-3. (58) 

The upper bound on macroscopicity in Section [HI B] 
then implies that the typical macroscopicity is necessar¬ 
ily small: Strongly geometrically entangled states cannot 
be macroscopic. 



FIG. 8: Macroscopicity of Haar-random states. Error-bars 
show one standard deviation, (a) Average macroscopicity and 
geometric entanglement, for N = 3,..., 23, no error bars 
shown, (b) Macroscopicity as a function of the number of 
qubits. Inset: Logarithmic plot with fit by an exponential de¬ 
cay. (c) Geometric entanglement. Blue solid line: Maximally 
possible value of geometric entanglement Eg = N—1. Dashed 
black line: Lower bound on geometric entanglement of ran¬ 
dom states, Eq. (58l. Error bars are not visible. Sample sizes 
(a-c) iV = 3 - 10 : 10®, iV = 11 - 15 : 10"^, iV = 16 - 23 : 10®. 
(d) Histogram for 10® Haar-random states of = 4 qubits, 
together with the upper bound of Section III B| (red solid line. 
Eg > 1) and the state |H) parametrized as in Fig.[^(red solid 
line and black dashed line, Eq < 1. 


3. Numerical results 

The expected behavior is reproduced by our numeri¬ 
cal data, shown in Fig. In agreement with the previ¬ 
ous argument, the geometric entanglement increases as 
a function of the number of qubits N (c), while the nor¬ 
malized macroscopicity decays (b). This decay is approx¬ 
imately exponential [logarithmic inset of Fig. |^b)], and 
we can safely state that even the unnormalized macro¬ 
scopicity M (Eq. 0) , which reflects the absolute size of 
the macroscopic component, decreases. The variances of 
both normalized macroscopicity and geometric entangle¬ 
ment decrease as well. The histogram in Fig.0d) shows 
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the distribution of states for = 4 in the {A4, £'G)-plaiie, 
together with the bounds in the regime Eq > 1 and the 
conjectured bounds in the realm Eq < 1. 

In conclusion, both the relative and the absolute size 
of the largest macroscopic superposition in Haar-random 
states decreases with N. As a consequence, Haar-random 
states are very geometrically entangled and feature little 
macroscopicity. 


C. Random symmetric states 


Colloquially speaking, Haar-random states are ex¬ 
tremely complex and do not allow any efficient descrip¬ 
tion . A state with large macroscopicity, on the other 
hand, can be approximated by a superposition of Eigen¬ 
states of the total spin operator [Eq. (52)], and thereby 
permits an efficient description. Hence, complexity and 
macroscopicity are mutually exclusive properties, and we 
cannot expect to encounter macroscopic superpositions 
in structureless ensembles. 

On the other hand, ensembles of random pure states 
that are less complex may feature higher values of 
macroscopicity. In particular, permutationally symmet¬ 
ric states constitute an ensemble of states with rather low 
geometric entanglement [5S1001 [Ml IM] : 


^^Gd^-sym)) <l0g2(iV + l), 


(59) 


due to the vastly reduced dimensionality A^ -|- 1 of the 
space of symmetric A^-qubit-states in contrast to the full 
Hilbert-space of size 2^. We choose the following ensem¬ 
ble of symmetric states: In the Dicke-state representation 
[Eq. (38)], the coefficients Cj are chosen to be normally 


distributed random variables (with normal real and imag¬ 
inary parts), and the resulting states are normalized. Fol¬ 
lowing this prescription, the ensemble is invariant under 
permutation-symmetric local unitary operations. 

The very different behavior of geometric entanglement 
for Haar-random and random symmetric states is evident 
comparing Figs. [^c) and l^c). Consistent with their low 
geometric entanglement, symmetric states feature excep¬ 
tionally high and robust macroscopicity. 

Expectation values of observables read 


N 


i'})\ 


0 dilTs) = ^ c*Cq{D^p\akiSia-i\Dj^ 
p,q=0 


(60) 


Since the Cp are chosen randomly and independently, only 
the summands with p = q will contribute to the average 
in the limit of many qubits N ^ oo. The only non-trivial 
expectation values that do not vanish on average are the 
two-qubit-correlations along the same axis, 




(61) 


Consequently, the matrix I4ym in Eq. (32) converges to 

N-1 


^ym — 


1 + 


(62) 


with obvious Eigenvalues, and, in the limit N —> oo, we 
therefore expect that the macroscopicity approaches 

Mdvk,))^^. (63) 

For finite N, off-diagonal non-vanishing correlations may 
contribute further to the fluctuations, which is why the 
average macroscopicity converges to l/-\/3 from above. 
This behavior is confirmed empirically in Fig. [^a,b), 
where the average value of macroscopicity for symmetric 
states is plotted against the average geometric entangle¬ 
ment (a) and the number of qubits (b). 


( 3 ) 
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FIG. 9: Macroscopicity of random symmetric states for a 
sample of 3000 random states. Error-bars show one stan¬ 
dard deviation, (a) Normalized macroscopicity against ge¬ 
ometric entanglement, for N = 2,..., 128. (b) Average 

normalized macroscopicity as a fnnetion of the nnmber of 
qubits N, the dashed black line shows the limiting valne 
l/-\/3. (c) Average geometric entanglement as a function of 
N, t he s olid line shows the theoretical maximum log 2 (N -|- 1) 
[Eq. ( |59| )] . (d) Average normalized largest and smallest Eigen¬ 
values A/(l -I- {N — l)/3) of the matrix I4ym for randomly 
chosen symmetric states. The largest Eigenvalue is directly 


related to the macroscopicity via Eq. (33l. The smallest 
Eigenvalue solves the minimization prob em that consists in 
finding the additive observable with the weakest fluctuations. 
For large N, the local spin orientation is rather irrelevant for 
experiencing large fluctuations, as long as all local spin mea¬ 
surements are performed along the same axis. 


Moreover, not only does the macroscopicity converge 
to a finite value, it is also very robust with respect to mis¬ 
alignment of spin-orientations: Since the smallest and 
largest Eigenvalues of I4ym [Eq. (^] converge to the 
same value 1 -I- (A^ — l)/3 [Fig. 0d) ], the spin-orient at ion 
becomes irrelevant in the limit Al —> 00 : Almost every 
additive observable for which the local spin orientations 
are all identical features macroscopic fluctuations on a 
random symmetric state. The equality of local spin ori¬ 
entations is crucial here: If these orientations were cho¬ 
sen randomly and independently, the expectation value 
would hardly fluctuate, since most Eigenvalues of the 
full variance-covariance-matrix V Eq. ( |I7| ) are typically 
small. 
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V. CONCLUSIONS AND OUTLOOK 

Many different approaches to entanglement eventually 
turn out to be equivalent, motivating the powerful con¬ 
cepts of entanglement monotone and entanglement mea¬ 
sure [55j . Our results emphasize that macroscopic entan¬ 
glement, as quantified by Eq. Q, should never be treated 
as a synonym for a measure of entanglement: In partic¬ 
ular, there are entangled states that feature vanishing 
macroscopic entanglement (Section 

Random physical states reflect this intricate relation¬ 
ship by their trajectory in the (Eg, Al)-plane (Section 


large geometric and small macroscopic entanglement. 
The typical size of macroscopic superpositions of ran¬ 
dom physical states grows, but not as fast as the system 
size - consequently, the normalized measure of macro- 
scopicity converges to 0 in the limit iV —>■ oo. Symmetric 
states are naturally much less geometrically entangled 
and much more macroscopic, and it remains to be stud¬ 
ied whether there are ensembles beyond symmetric states 
for which the actual spin orientations in the definition of 
the additive observables [Eq. Q] are irrelevant. Such 
ensembles would be experimentally valuable due to their 
robustness. 

Further quantitative insight in the relation between 
macroscopic and geometric entanglement is desirable. A 
general bound on geometric entanglement as a function 
of macroscopicity (or vice versa) seems hard to obtain, 
since both quantities are defined via a maximization pro¬ 
cedure. We believe nevertheless that our bound in Sec¬ 
tion jlTT^ can be improved considerably, and that a proof 
for the extremality of |S(0, e)) can be found. We did 
not find any relationship between the closest separable 
state |(()i, ((> 2 ,..., (^at) and the maximizing spin orienta¬ 
tion {di,..., aAf}; such deeper connection would be valu¬ 
able. 


IV A), converging to Haar-random states, which feature 


Control schemes that optimize multipartite entangle¬ 
ment implicitly exploit the typicality of entangled states 
within the ensemble of pure states [SB]. Our results sug¬ 
gest that control strategies that aim at a macroscopically 
entangled target state will not only be affected by deco¬ 
herence, but the unitary evolution also needs to be tai¬ 
lored in a much more precise way: While the manifold of 
states that feature high geometric entanglement is very 
large, this is not true for macroscopic states. 

Finally coming back to our proposed analogy [Fig. [^, 
our results suggest that macroscopically entangled states 
play the role of four-leaf clover: They do not appear spon¬ 
taneously after some random process, but only as the 
result of some meticulously designed artificial evolution, 
such as in a quantum computer m- Further investiga¬ 
tions of other ensembles of pure quantum states, such 
as canonical thermal pure states [58] and random matrix 
product states [591160] , will eventually rigorously confirm 
or dismiss the analogy. 


Acknowledgements 


C.-Y.P., M.K., and H.J. were supported by the Na¬ 
tional Research Foundation of Korea (NRF) grant funded 
by the Korea Government (MSIP) (No. 2010-0018295). 
M.C.T and K.M. acknowledge financial support by the 
Danish Council for Independent Research and the Villum 
Foundation. M.C.T. was financially supported by the bi¬ 
lateral DAAD-NRF scientist exchange programme. C.- 
Y.P. thanks Jinhyoung Lee and Seokwon Yoo for granting 
access to the Alice cluster system at the Quantum Infor¬ 
mation Group, Hanyang University. The authors thank 
Christian Kraglund Andersen, Ralf Blattmann, Fliska 
Greplova, Qing Xu and Jinglei Zhang for valuable com¬ 
ments on the manuscript. 


[1] E. Schrodinger, in Proc. Cambridge Phil. Soc., p. 555, 
1935. 

[2] E. Schrodinger, Naturwissenschaften 23, 823 (1935). 

[3] F. Frowis and W. Diir, New J. Phys. 14, 093039 (2012). 

[4] A. Shimony, Ann. N. Y. Acad. Sci. 755, 675 (1995). 

[5] H. J. Briegel and R. Raussendorf, Phys. Rev. Lett. 86, 
910 (2001). 

[6] A. Leggett, Prog. Theor. Phys. Suppl. 69, 80 (1980). 

[7] A. Leggett, J. Phys. Condens. Matter 14, R415 (2002). 

[8] Y. Chen, J. Phys. B At. Mol. Opt. Phys. 46, 104001 
(2013). 

[9] C. Myatt, B. King, Q. Turchette, C. Sackett, D. Kielpin- 
ski, W. Itano, C. Monroe, and D. Wineland, Nature 403, 
269 (2000). 

[10] W. H. Zurek, Rev. Mod. Phys. 75, 715 (2003). 

[11] M. Schlosshauer, Rev. Mod. Phys. 76, 1267 (2005). 

[12] A. Buchleitner, M. Tiersch, and C. Viviescas, Eds., En¬ 
tanglement and Decoherence: Eoundations and Modern 
Trends Lect. Notes Phys. 768 (Springer, Berlin, Heidel¬ 


berg, 2009). 

[13] H. Brener and F. Petruccione, The Theory of Open Quan¬ 
tum Systems (Oxford University Press, Oxford, 2006). 

[14] S. Raeisi, P. Sekatski, and C. Simon, Phys. Rev. Lett. 
107, 250401 (2011). 

[15] P. Sekatski, N. Sangouard, and N. Gisin, Phys. Rev. A 
89, 012116 (2014). 

[16] P. Sekatski, N. Gisin, and N. Sangouard, Phys. Rev. Lett. 
113, 090403 (2014). 

[17] K. Huang, Statistical Mechanics (Wiley, Hoboken, NJ, 
1987). 

[18] T. Morimae, Phys. Rev. A 81, 010101 (2010). 

[19] A. Streltsov, H. Kampermann, and D. Brufi, Phys. Rev. 
A 84, 022323 (2011). 

[20] D. Gross, S. Flammia, and J. Eisert, Phys. Rev. Lett. 
102, 190501 (2009). 

[21] M. Bremner, C. Mora, and A. Winter, Phys. Rev. Lett. 
102, 190502 (2009). 

[22] W. Diir, G. Simon, and J. Cirac, Phys. Rev. Lett. 89, 









13 


210402 (2002). 

[23] G. Bjork and P. G. L. Mana, J. Opt. B Quantum Semi- 
classical Opt. 6, 429 (2004). 

[24] A. Shimizu and T. Morimae, Phys. Rev. Lett. 95, 090401 
(2005). 

[25] J. I. Korsbakken, K. B. Whaley, J. Dubois, and J. I. 
Girac, Phys. Rev. A 75, 042106 (2007). 

[26] G.-W. Lee and H. Jeong, Phys. Rev. Lett. 106, 220401 

( 2011 ). 

[27] Y. Y. Xu, Phys. Scr. 88, 055005 (2013). 

[28] H. Jeong, M. Kang, and H. Kwon, Opt. Commun. 337, 
12 (2014), 1407.0126. 

[29] T. Farrow and V. Vedral, Opt. Commun. 337, 22 (2015). 

[30] S. Nimmrichter and K. Hornberger, Phys. Rev. Lett. 
110, 160403 (2013). 

[31] A. Laghaout, J. S. Neergaard-Nielsen, and U. L. Ander¬ 
sen, Opt. Commun. 337, 96 (2014), 1405.5062. 

[32] F. Frowis, N. Sangouard, and N. Gisin, Opt. Commun. 
337, 2 (2015). 

[33] D. M. Greenlrerger, M. A. Horne, and A. Zeilinger, Go¬ 
ing beyond bell’s theorem, in Bell’s Theorem, Quan¬ 
tum Theory and Conceptions of the Universe, edited by 
M. Kafatos, Kluwer, 1989. 

[34] A. Shimizu and T. Miyadera, Phys. Rev. Lett. 89, 270403 
(2002), 0203106. 

[35] T. Morimae, A. Sugita, and A. Shimizu, Phys. Rev. A 
71, 032317 (2005). 

[36] D. Baguette, T. Bastin, and J. Martin, Phys. Rev. A 90, 
032314 (2014). 

[37] H. Barnum and N. Linden, J. Phys. A Math. Theor. 34, 
6787 (2001), 0103155. 

[38] J. Eisert and H. Briegel, Phys. Rev. A 64, 022306 (2001). 

[39] H. Mine and M. Marcus, Permanents (Cambridge Uni¬ 
versity Press, Cambridge, 1984). 

[40] J. Martin, O. Giraud, P. A. Braun, D. Braun, and 
T. Bastin, Phys. Rev. A - At. Mol. Opt. Phys. 81, 062347 
( 2010 ). 

[41] R. Hiibener, M. Kleinmann, T. C. Wei, C. Gonzalez- 


Guillen, and O. Giihne, Phys. Rev. A 80, 032324 (2009). 

[42] T. J. Volkoff and K. B. Whaley, Phys. Rev. A 90, 062122 
(2014). 

[43] A. Hamma, S. Santra, and P. Zanardi, Phys. Rev. A 86, 
052324 (2012). 

[44] A. Hamma, S. Santra, and P. Zanardi, Phys. Rev. Lett. 
109, 040502 (2012). 

[45] W. Wootters, Found. Phys. 20, 1365 (1990). 

[46] A. Sugita and A. Shimizu, J. Phys. Soc. Japan 74, 1883 
(2005). 

[47] M. C. Tichy, F. Mintert, and A. Buchleitner, J. Phys. B 
At. Mol. Opt. Phys. 44, 192001 (2011). 

]48] K. Zyezkowski and M. Kus, J. Phys. A. Math. Gen. 27, 
4235 (1994). 

[49] M. Tiersch, Benchmarks and Statistics of Entanglement 
Dynamics, PhD thesis. University of Freiburg, 2009. 

[50] M. Ledoux, The Concentration of Measure Phenomenon, 
Mathematical Surveys and Monographs Vol. 89 (Am. 
Math. Soc., 2001). 

[51] N. Ullah, Nucl. Phys. 58, 65 (1964). 

[52] H. Venzl, A. Daley, F. Mintert, and A. Buchleitner, Phys. 
Rev. E 79, 056223 (2009). 

[53] D. J. H. Markham, Phys. Rev. A 83, 042332 (2011). 

[54] M. Aulbach, D. Markham, and M. Murao, New J. Phys. 
12, 073025 (2010). 

[55] M. B. Plenio and S. Virmani, Quant. Inf. Comput. 7, 1 
(2007). 

[56] F. Platzer, F. Mintert, and A. Buchleitner, Phys. Rev. 
Lett. 105, 020501 (2010). 

[57] A. Shimizu, Y. Matsuzaki, and A. Ukena, J. Phys. Soc. 
Japan 82, 054801 (2013). 

[58] S. Sugiura and A. Shimizu, Phys. Rev. Lett. 108, 240401 

( 2012 ). 

[59] S. Garnerone, T. R. de Oliveira, and P. Zanardi, Phys. 
Rev. A 81, 032336 (2010). 

[60] S. Garnerone, T. R. de Oliveira, S. Haas, and P. Zanardi, 
Phys. Rev. A 82, 052312 (2010). 



